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Abstract 

The consecutive microscopic solution is presented of the problem of tunneling 
of a particle through a potential barrier. The method is applied to the a- and 
proton decay of compound systems formed in fusion reaction. Appearance of the 
peaks in the spectrum of emitted particles is predicted. The peaks correspond 
to quasistationary states inside the potential barrier. 



1 Introduction 



Prefission emission of alphas and light charged particles provides with a source of 
information about the time scale in the fusion-fission reaction. On the other side, 
theoretical description of the processes is usually made in terms of the inverse cross- 
section. Such an approach assuming time reversibility of the process leaves out of the 
scope a possibility of its profound experimental check. This is in contrast with a number 
of indications that violations of the reversibility may arise due to back-transparency of 
the inner slope of the potential barrier in the ingoing channel [1], or different response of 
the nuclear surface on the interaction with the emitted and the same incoming particle, 
or due to temperature effects on the barrier distribution [2], as in the ingoing channel 
experimental fit of the optical model parameters is only possible for cold nuclei [3]. 

Moreover, traditional decay theory deals with tunnelling through a barrier of a 
particle which is in the quasistationary state. This does not involve important cases 
when a virtual proton or cluster is between the quasistationary states, as in alpha decay 
from compound systems formed in fusion reaction. 

Our approach allows one to calculate the decay width at any energy of the emitted 
particle. Strong resonance effects are, specifically, predicted in alpha spectra from 
compound systems produced in heavy-ion collisions. 

2 Formalism 
2.1 a decay widths. 

Let \i/(ri, . . . , r^) be a wave function of the source nucleus with the mass number A. 
First, it can be expressed in terms of the channel wave-function basis, as products 
of the wave-function of the daughter nucleus (y9„(ri, . . . , rA-4) and the a-particle w.f. 
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where we selected the angular momentum L of the relative motion of the a particle 
in the nucleus, rfl may be treated as a wave-function of the relative motion of the 
a cluster in the mother nucleus which evidently turns out to depend on the relative 
coordinate R. In simple cases of pure configuration the expansion coefficients CnkL are 
reduced to genealogical coefficients. 

Let then the nucleus make a transition i f. As a result, in the exit channel we 
observe the system in a state which is described by a wavefunction as a superposition 
of the plane wave and ingoing spherical wave |1] at large a-nucleus distances R: 

V'/p(ri, . . . ,rA) ~ VP/(ri, . . . ,r^_4)5(p(rA_3, . . . ,rA) , (2) 
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In eq. (3), Qp is the channel wave-function, which is the eigen function of the a-nucleus 
Hamiltonian with an appropriate mean- field single particle potential Ua{R)'- 
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Bp = pV2M„ . (5) 

Furthermore, taking into account the asymptotics (3), the wavefunction jFp(R) can be 
expressed in terms of the spherical harmonics in a usual way: 

oo 
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To find a transition amplitude, one has to change to the coordinate system of the 
exit channel |/p). The transformation of eq. (1) then conventionally reads as 

p 

the re-expansion coefficients giving the transition amplitude under consideration. This 
way is similar to that found by Migdal when solving his classical problem of shake of 
an atom in /? decay jl]. Substituting eqs. (1) together with (2), (3) and (6) into eq. 
(7), we arrive at the following expression for the transition amplitude: 
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Taking into account that the wave functions JFp are normalized at 1 particle in a unit 
volume, with the flux v = p/M^, we obtain from (8) the following expression for the 
decay probability per a unit time: 

rp-^ = |M/pP- (9) 

Inserting Mfp from eq. (8) into eq. (9) and integrating over all the angles of emission 
within An, we arrive at the following final expression for the decay width: 

r« - ^ = 4vrM>5^ f; icM'mrn' . (lo) 
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3 Method of numerical solution. Eigenvalues 



a-nucleus potential is characterized by a Coulomb barrier, which is high enough, to 
form quasibound states inside the barrier (Fig. H]). 

These would be usual eigenstates, if the barrier were infinitely broad. The values, 
however, go over the resonances on the continuum background, whenever the pene- 
trability of the barrier is taken into account. Coupling to the continuum causes the 
energy shift and broadening of the eigenstates. Affected eigenvalues can be determined 
as follows. 

The Schrodinger equation for an a-particle in the field of a nucleus reads as follows: 
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and the Coulomb potential was taken into account as due to the sharp-edge charge 
distribution: 
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with Rq being the nuclear radius. 
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Figure 1: a-nucleus potential for the system of ^'^^La + a with the zero angular momentum. 



On the radius segment between the origin r = and the first turning point Rd'- 
< r < Rci, eq. (1) was integrated numerically with the initial condition 

^(r) ^ . (15) 

r ^ 

General solution of the Schrodinger equation under the barrier is a linear combina- 
tion of the two linearly independent solutions. One of them exponentially vanishes, and 
the other exponentially increases with increasing R. The coefficients can be obtained 
by sewing the functions at the internal turning point. In principle, the eigenvalues 
may be obtained from a condition that the coefficient for the exponentially increasing 
solution vanishes. In the first approximation, this can be achieved by sewing to the Airi 
function j4j. Actually, the eigen solutions were obtained by numerical integration from 
the external turning point i?c2 towards the internal one, with somewhat an arbitrary 
initial condition 

yiR,2) = l, y'{Rc2) = -0.25. (16) 

The derivative in eq. (6) is negative, as the solution is assumed to exponentially de- 
crease under the barrier. In the course of integration, only the right solution survives, 
which exponentially increases with decreasing R under the barrier, the other exponen- 
tially vanishes, in so far that the eigenvalue obtained practically very weakly depends 
on the concrete numbers in eq. (jlfjj) . 

In general case, the fundamental set was obtained by numerical integration from 
i?c2 to Rci with two different initial conditions: 

y{Rc2) = l, y'{Rc2) = ±l. (17) 

After sewing at r = R^i, the resulting solution increases under the barrier (see below 
Figs. El and 0} if not an eigenstate, in contrast with the behavior of each of the fun- 
damental solutions. This demonstrates mathematical correctness of the method. For 
the numerical integration, the Runge-Kutta-Nystrom method was used. The Shtermer 
method was also tried, with essentially the same results. Behind the barrier, the both 
solutions oscillate. 

4 Results and discussion 

Calculations were performed with the Saxon- Woods potential ()1H|1 . with the parameters 
Vq = 100 MeV, s = 2.3 Fm, c = 1.2A^^^ Fm. Representative wavefunctions for various 
energies are presented in Figs. |2Hllfor the system a + ^'^^La, L = 0. 

Fig. 121 answers the eigenvalue of = 10.79 MeV. Corresponding wavefunction 
has a large amplitude inside the barrier. Therefore, the overlapping integral is also 
expected to be large in this case. 

In Figs. El and m we present the wavefunctions aside the resonance, for the energies 
of 11 and 14 MeV, respectively. The wavefunctions are normalised at S{p — p'). These 
figures are in drastic contrast with the resonance one, presented in Fig. |21 The ampli- 
tude of the wavefunction within the barrier is much smaller than outside. As a result, 
the overlapping integral is expected to be small in the nonresonance case, depressing 
nonresonance a decay. 
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Figure 2: a wavefunction for the system of ^■^^La + a with the a energy of 10.79 MeV. 

Finally, in Fig. Elwe present the calculated line T{E). That has a typical resonance 
shape with a half-width of around 200 eV. Therefore, spectrum of the subthreshold a 
particles turns out to be modulated, directly indicating the resonance states inside the 
barrier. Confirmation of this effect in experiment would really mean discovering new 
physics. 

In heavy-ion collisions, this effect may be smoothed by mixed multipolarities. The 
effect must also manifest itself in usual a- or proton decay, specifically, of nuclei far from 
the drip line. In this case, set of the allowed L values is usually not large. Moreover, 
a partial wave with a certain L may make predominant contribution, which can be 




Figure 3: a wavefunction for the system of ^'^^La + a, with the a energy of 11 MeV. 
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Figure 4: a wavefunction for the system of ^■^^La + a, with the a energy of 14 MeV. 

exploited for direct check of the theory presented herein. This study is planned to be 
made separately, in due course. 
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Figure 5: Profile of the a decay line in an anticipated sub-barrier spectrum, half-width being around 200 
eV. (The system is ^■^^La + a, for the alphas emitted with L — 0.) 
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